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Classical escape in 2D Hamiltonian systems with the mixed state has been studied numerically and 
analytically. The wide class of potentials with the mixed state is presented by polinomial potentials. 
In potentials, where the mixed state could be realized, i.e. the phase space contains regions of both 
regular and chaotic motion, escape problem has a number of new features. In particular, some local 
minima become a trap with number of particles depending on energy and other values that characterize 
the ensemble of particles. Choosing the form of initial ensemble one chooses the set of parameters that 
determine the number of trapped particles. 

1 Introduction 

Escape of the trajectories from the locahzed regions of phase or configuration space is 
an important topic in dynamics and describes the decay of metastable states in many 
areas of physics, as for instance chemical and nuclear reactions, atomic ionization and 
others. The problem has a rich history and a number of realizations in different systems. 
Almost a century ago Sabine |lL| had considered the decay of sound in concert halls and 
later Legrand and Sornette had shown that problem is equivalent to the escape one. 
Corresponding decay rate is / a{s)ds and a{s) is the absorption coefficient at coordinate 
s of billiard boundary, a{s) = 1 at the opening of width A and a{s) = elsewhere. 

Another application of the escape problem links to the nondestructive monitoring of 
the system [3j . If some system is connected to the surroundings only via small opening in 
its boundary, it became possible to understand the dynamics of the system by exploring 
the escaping particles. So the natural question arises: how does escape law depend on 
the character of the motion? For strongly chaotic systems exponential decay is expected 
[H El E]. Bauer and Bertsch [1] considered the escape of the particles through the small 
opening in the billiards boundary. When exploring the regular billiard, i.e. rectangular 
without the scattering center, power law emerges at long time. Qualitative understanding 
of the mechanism of power tails generation is given in [7j. 

For rectangular billiard with circular scattering center inside the decay of the initial 
ensemble of N{0) particles is exponential and by simple considerations one could obtain 
the corresponding decay rate: 
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N(t) = N(0) exp(-at), a = (1) 

TcAc 

here p is the particles momentum, A - the width of the opening and Ac - the area of 
the biUiard. As it would be shown further, exponential decay is a common feature of the 
purely chaotic systems. 



2 The mixed state 

Passing from billiards to the potential systems broadens the number of possible situations. 
One-well potential is the simplest case for considering the escape. Zhao and Du [8J have 
explored the escape from Henon - Heiles potential: 

UHH{x,y) = ^^ + xy'-^ (2) 

At the energies E > 1/6 trajectory could leave the potential well through one of the 
three openings, placed symmetrically. Numerical simulation, performed by Zhao and Du, 
has shown that escape follows the exponential law. At over-saddle energies phase space 
of the Henon - Heiles Hamiltonian is almost homogeneous and motion is chaotic. Using 
this, the escape rate could be derived and it fits numerical results with high accuracy. 
This situation is similar to the escape from chaotic billiards. Another study of the escape 
from Henon - Heiles potential was performed by [9]. 

In contrast to billiards and Henon - Heiles potential, some potentials have highly 
inhomogeneous phase space, that consists of macroscopically significant components of 
both regular and chaotic motions. One wide class of systems with inhomogeneous phase 
space is represented by multi-well potentials. We will focus our research on the escape in 
such potentials. The preliminary results are presented in [TT]. Regularity-chaos transition 
in multi-well potentials has a distinctive feature, which consists in the difference of critical 
energies in different local minima. This leads to the different (either regular or chaotic) 
regimes of motion in different local minima at the same energy, i.e. the ratio of chaotic 
trajectories in some local minimum significantly differs from the ratio in other minima. 
Such kind of the dynamics is called the mixed state . It is important to mention, that 
critical energies lie below the saddle energy, i.e. the energy above which local minima are 
no more separated. 

We will demonstrate the mixed state in two representative examples of 2D multi-well 
potentials: the lower umbilic catastrophe D^, 

Ud, = 2ay^ - + xy^ + ^x^ (3) 
for a = 1.1 and the quadrupole nuclear oscillations potential (QO), 

Uqo{x, y, W) = ^^^^ + V - + ix' + yr (4) 

for W = 18. Z^5 potential has two local minima and three saddles and it is the simplest 
potential, where where mixed state is observed. Fig. [1] shows the Poincare sections for dif- 
ferent energies in the considered potentials. It demonstrates the evolution of dynamics in 
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different local minima. At low energies motion has well marked quasiperiodic character in 
both minima. As energy grows, gradual regularity-chaos transition is observed. However, 
changes in features of trajectories, localized in different minima, are sharply distinct. In 








-0.1 0.1 0.2 0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 

X X X X 

Figure 1: Level lines and Poincare sections for D5 (upper row) and QO potentials at different energies 

the left minimum significant fraction of trajectories becomes chaotic already at about a 
half of the saddle energy, and at near saddle energy almost all initial conditions result 
in chaotic trajectories. In the right minimum at the same energy the motion remains 
regular and this situation is preserved up to the saddle energy (we will call this minimum 
the regular one for simplicity). Moreover, at energies above the saddle one the phase 
space is still divided on chaotic and regular components, but they are not separated in 
configuration space. 

Earlier we have shown that the mixed state opens new possibilities for investigations of 
quantum manifestation of classical stochasticity [12] • Aim of the present work is to study 
the classical escape from separated local minima, realizing the mixed state. We show that 
escape from such local minima has all above mentioned properties of the decay of chaotic 
systems and also a diversity of principally new features, representing an interesting topic 
for conceptual understanding of chaotic dynamics and for applications as well. We are 
interested in both the "first passage" effects and dynamical equilibrium setup for the finite 
motion (for example, in QO potential). It is important to stress, that though we study the 
process of escape from concrete local minimum, the over-barrier case of the mixed state 
has specific memory: general phase space structure at supersaddle energies is determined 
by the characteristics of the motion in all other local minima. 
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Figure 2: Poincare section for D5 potential &t E — 1.0. Initial conditions for quantum computations 
are presented 

3 Decay of the uniformly distributed ensemble 

At the energies above the saddle, i.e. E > Es, different components are not separated 
in the configuration space. Fig. [2] represents the Poincare section for D5 potential at 
supersaddle energy. The "chaotic sea" stretches on whole accessible area, while regular 
island in the right well is localized. This means that, been initially localized in right well, 
chaotic trajectories could leave the well and regular ones remain trapped, i.e. the decay 
of the mixed state occurs. Therefore we will explore the escape of the particles from the 
right well of and peripheral wells of QO potential. 

Numerical simulation of the escape process in these potentials implies three steps. At 
the first stage we select initial distribution of the particles inside the well. Then direct 
numerical integration of the equations of motion for all particles is performed and we 
extract N{t)/N{0) — relative number of particles in the well. Using this function one 
can calculate escape rate and part of the trapped trajectories. One remark should be 
made about the first step of the numerical simulation. Initial distribution, in general, 
determines the ratio between regular and chaotic trajectories and hence it should be 
physically motivated. One chance is to distribute particles uniformly in all classically 
allowed configuration space and another — to put all particles at the same point. Second 
case emulates injection of particles to the well. In both cases momentum is calculated 
using the energy conservation (in the second it will be the same for all particles) and 
its direction is uniformly distributed in [0,2-7?]. These initial distributions present quite 
simple extreme cases of real distributions. Uniform distribution will be illustrated with 




Figure 3: Decay law N{t)/N{0) for D5 potential at different energies. E, = OA + 0.2k, k = O.A 

Ud^ and "point" distribution - with Uqq- 

Phase space density for uniform initial distribution is 

where S{E) is the area of classically allowed space: 



S{E)= J dxdye{E -U{x,y)) (6) 

X>Xg 

Numerical simulation reveals some substantial features of escape process for this initial 
distribution. Fig. [3] demonstrates the normalized number of particles in the well as a 
function of time. Decay law has three important features: 

• saturation at t ^ 00: 

N{t ^ 00) = pooA^o (7) 

Because of uniform initial distribution, quantity poo is a relative phase volume, 
occupied by trapped trajectories. 

• initial linear decrease - from to some t{E): 

N{t)/No = 1 - a'-^h (8) 
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Figure 4: Poincare section for trapped trajectories. These trajectories form regular island, i.e. they are 
regular 

• exponential decrease at t > t{E): 

N{t)/No = Poo + Cexp{-a^'h) (9) 

Fig. m presents the Poincare section for trapped trajectories. Obviously, regular island 
in the right well is formed by trapped, i.e regular, trajectories. Function Poo{E) is demon- 
strated on Fig. O This is decreasing function, as expected from physical considerations — 
number of regular trajectories decreases as energy increases. 

It is interesting to mention, that there exists correlation between poo{E) and relative 
area of regular island in Poincare section p. While relative part of trapped trajectories 
is linked to the volume of four-dimensional phase space, occupied by regular trajectories, 
Poincare section is two-dimensional and thus there is no argument to expect the precise 
coincidence of p and poo{E). 

To calculate the relative area of regular island in Poincare section we first determine 
the border of the island through numerical integration of equation of motion and calculate 
the area inside it. Then this area was divided by total area, determined by the conditions 
X > 0, > 0. In spite of topological inequivalence, p and Poo{E) are very close to each 
other. This means that one could determine and control the part of trapped trajectories 
using only the Poincare section. 

Linear part of the decay law is more pronounced in comparison with pure ensemble 
[8]. In the time interval between and t{E) decay has the form ([8]). Thus, there are two 
quantities describing linear decay — its duration t{E) and corresponding escape rate a^'^. 

After linear part, at t > t{E), decay law has exponential form (Q. It's important that 
(E]) is not a linear approximation of ([9]), and linear decay has independent nature. This 
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Figure 5: p, the relative area of the regular island in Poincare section, is shown by squares and Poo{E) 
is shown by circles 



means that p^o + C 1. Moreover, decay law of the form ([H]) precisely works up to time 
t ~ t{E) and in this interval it differs substantially from corresponding exponential law 
Poo + (1 - Poo) exp(-a(^)t). 

Let's now calculate t{E). From analysis of numerical calculation one can derive that 
this time corresponds to the time of one-dimensional, along y = 0, motion from the saddle 
to opposite side of the well and back. Thus t{E) has the form: 



rD,{E) 



TQo{E) = 12 



dx ^ j<[ 

\v\ ~ E-i 
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(10) 



where K{k) is full elliptic integral of the first type and Es = 1/12^ — saddle energy in QO 
potential at = 18. 

Such nature of t{E) and analysis of linearly escaping trajectories allow to conclude 
that linear decay corresponds to the escape of trajectories that move along y = and 
cross the well not more than two times. Correspondingly, trajectories which initially move 
toward saddle along y = escape first, and then escape occur for trajectories moving to 
the opposite part of the well. 

a'-'^ could be calculated via averaging of the flow through the saddle: 

r7r/2 
l-n/2 



a^^\E) = p{E) ( dy d9 \v\ cos 6 

J J-it/2 



(12) 



x=xs 
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This procedure is the same as in [8]. Using the density ([5]) and integrating the (fT2l) we 
obtain the expression for hnear escape rate: 




Exponential decrease of N{t)/N{0), as it could be understood from the analysis of escaped 
trajectories, corresponds to the leaving of sticking orbits, i.e. those chaotic trajectories 
which moved in the vicinity of the regular island in the Poincare section. 

Energy is the parameter which determines the part of trapped trajectories for the 
uniform ensemble. Changing the energy of ensemble one could trap the given number of 
particles. 
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Figure 6: Initial ensemble for extraction of asymptotic distribution. 

To illustrate the splitting of initial ensemble on regular and chaotic components one 
needs to plot the ensemble in {x,y) plane at time t ^ t{E). It is convenient to plot the 
asymptotic evolution of the initial ensemble in QO potential because of finite character 
of motion in it. At t = particles are equally distributed between peripheral minima. 
Fig. [6] represents this initial ensemble. 

The ensemble splits during the evolution in time — regular trajectories remain trapped 
in peripheral minima, while chaotic trajectories cover almost entire accessible configura- 
tion space. Integrating equations of motion for all trajectories in the initial ensemble for 
enough long time (in fact — for time much greater than the typical escape time) we ob- 
tain the particles positions corresponding to asymptotic distribution. We have calculated 
the asymptotic distribution for energy E = l.5Es {Es = 1/12"^ - saddle energy in QO 
potential), while integration was performed for t = 150 (for this energy t{E) = 28.395). 
Fig. [7| shows the asymptotic distribution of trapped particles. At enough large time these 
particles tend to accumulate closer to the center of the well. 
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Figure 7: Asymptotic distribution of trapped particles. 
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Figure 8: Asymptotic distribution of free particles. 



The corresponding distribution of free particles is represented at fig. [HI As it was 
mentioned above, free particles cover entire central minimum and, according to their 
character, the area near the plane y = in peripheral minima. If free particles are 
removed in some way after their leave the peripheral minima (this is, of course, correct 
for right well of and any well of same topology), we obtain pure regular ensemble 
inside the well, i.e. the initial mixed ensemble splits. 



4 Point ensemble 

For simplicity we will consider the point distribution with y = 0. Thus the governing 
parameter is xq — x-coordinate of the point of injection. Corresponding density has the 
form: 

S {x - xo} S {y} 6 {p - ^2iE - U{xo, 0))} 

p [x, y, p, (f) = / / . = (1^) 

27ry^2(E-f/(xo,0)) 
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Numerical procedure is identical to that for uniform ensemble, but for distribution ( fTSl) 
one could made general conclusion about the character of decay even without numerical 
integration. 

For this one needs to consider the representation of initial distribution on the Poincare 
section: 



^ ^ 5 {x - xo} e {p, - J2{E - Uixo, 0))} {J2iE - Uixo, 0)) + , , 

pps{x,Px) = ^ , / = (16) 

2J2{E-U{xo,0)) 



where 9 {a} is a step function. Lets denote: 



^(max)^^2(E-f/(Xo,0)). (17) 

Thus in the Poincare section initial ensemble occupies the interval x = Xo,Px E — p^™^^-* , p^™^^-* 
Edges of this interval correspond to the momentum directions 99 = vr, 0. In the over-barrier 
case this interval crosses the regular island in the points {xo,pi^^^^) and (xq, —p^^^^)- Ob- 
could not leave the well. Value 



ni'^eg) ^(reg) 
fx 1 fx 



viously, particles with p^ in the interval 

V?max(-E,Xo) = 2arccos 



(max) 
kPx 



defines the cone of directions which could leave the well. In other words, particles with 
Py which is greater than some maximum value are trapped. Thus first conclusion about 
decay of point ensemble implies the existence of escape cone and this feature reveals the 
role of transversal momenta in the escape process. 

The second feature consists in the fact that decay begins at the time ri, which corre- 
sponds to the time of motion of particle with momentum px = —p from point Xq to the 
saddle: 

xo xo 



n= —= , (19) 

I Vx I ^2{E-Uix,0)) 

Moreover, escape is a two-stage process due to existence of escape cone. At the second 
stage the particles, moving toward the well boundary in the escape cone, leave. 

Numerical procedure implies the determination of the ip^axi^o, E) , ti{xo,E), and 
N{t = 00), i.e. the number of particles in the well. Due to uniform distribution of 
momenta directions one has the relation: 

Poo(Xo, E) = = 1 (20) 

We will illustrate the above consideration for QO potential. Energy is normalized to 
the saddle energy, Es = 1/144^: E = sEs- Initial ensemble is localized in peripheral 
minimum with Xmin = 1/6 and corresponding saddle xs = 1/12. Fig. represents the 
normalized number of particles in the well as a function of time for different values of 
energy and xq = 0.16. 



10 



N{t]/N{0) 




60 



80 



100 



Figure 9: Normalized particle number in the well for point ensemble at different energy and xq = 0.16 



Numerically obtained ti could be compared to analytical value: 



Xo 



1/12 



dx 

Px 



xo 



dx 



1/12 



2{eEs-UQo{x, 0)) 



(21) 



Numerical and analytical value of first escape time are presented at fig. [TO] 

The most interesting question is the correspondence between escape cone angle, num- 
ber of trapped particles and linear part of regular island in Poincare section. FigJTTl 
represents the quantities 



Pps 



(max) / \ 

Px [x,e) 



(22) 



and 



= l-cos(V'max(a;,e)/2) (23) 

First is a linear part of regular island in the section and second is the corresponding 
expression through the escape cone angle. 

This angle could be determined during numerical simulation. Thus, Poincare section 
could be used to determine the angle of the escape cone. On the other hand, this angle 
is connected with the part of trapped trajectories. To demonstrate this connection one 
needs to compare 

</2max(a;, e) 



1 



(24) 



with Poo- Corresponding data are presented at fig. [121 This analysis allows to conclude 
that point ensemble differs substantially from the uniform when considering the escape. 

In uniform ensemble Poincare section gives only estimate (nevertheless very accurate) 
of trapped particles number, while in the point ensemble one could calculate not only part 
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Figure 10: 
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Figure 12: and p^p for different energies and injection points 



of trapped particles, but also the escape cone angle, using only Poincare section. Point 
ensemble allows dual control of trapped particles part using the energy and injection 
point. One even needs not to compute the entire Poincare section, but only the boundary 
of regular island, and then use the relation: 



The procedure of calculating the part of trapped particles, thus, implies three steps. First 
of all one needs to calculate the boundary of the regular island in the Poincare section. 
After that, the •p^^'^^^ and p^™^^ should be determined. Now value poo could be calculated. 

5 Quantum escape problem 

Now we consider the over-barrier decay of the mixed state from the quantum-mechanical 
point of view. In semiclassical limit, temporal evolution of an initial state in form of the 
minimum uncertainty Gaussian wave packet 



represents the quantum analogue for classical motion of a point particle with initial con- 
dition (xo, yojPxcPi/o)- The quantum-classical correspondence between the wave packet 
motion and the classical trajectory preserves quite a long time until the wave packet 
spreads. 



Poo 



2 arccos 




(25) 



TT 
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We made numerical simulations for time evolution of the Gaussian wave packets for 
three different initial conditions in the potential (3) with a = 1.1. The following wave 
packets parameters 



h = 0.01, ar 



0.07, a,, 



\ 



0.06, xo = V2, 



yo 



and Pq = '^Pxo + P^o = 2 were the same for all the three initial states, so all of them 
started from the right local minimum of the potential with initial energy E = 1, which is 
twice higher than the barrier height. (Recall that in the potential (3) Emin = — 1 and 
Eg = 0). The only difference was in the direction of initial momentum: we considered 
the cases 



TT TT 

«• 4- 2- 



The initial condition with (fQ = it/ 2 falls near the center of the principal stability 
island and therefore it corresponds to regular classical trajectory trapped in the right 
potential well. Quantum autocorrelation function 



P{t) = j dxdy^{t = O)^(t) 



(26) 



manifests quasiperiodic nature of the corresponding classical trajectory: the sharp peaks 
indicate periodically repeated recurrences of the wave packet to the initial state, and rather 
high amplitude of the peaks shows that there is almost no spreading of the wave packet. 
The reason of that slow spreading is the fact that even for highly over-barrier energies the 
regular classical motion near the potential minimum is still pretty close to that in two- 
dimensional harmonic potential, as it is seen for example in the characteristic structure 
of the stability island on the Poincare section. For the considered case of the right local 
minimum in the potential the corresponding quadratic (harmonic) approximation 
reads 



U{x,y) 



ujlix 



- V2f + u:y 



(27) 



with ujx = 2 and Uy = 2 y a + 1/ ~ 2.7 for a = 1.1. Therefore the trapped quasiperiodic 
trajectories actually are very close to the Lissajoux figures — simple superpositions of 
harmonic oscillations in perpendicular directions. Remarkably, the considered gaussian 
wave packet coincides with the coherent state for the harmonic oscillator potential flTTl) — 
the exact non-spreading solution of time depending Schrodinger equation — that is why 
it is "almost" non-spreading near the minimum of the potential. As the considered 
initial condition (fo = 7t/2 leads to almost one- dimensional motion (along the y-axis), only 
one frequency Uy ~ 2.7 is manifested in the autocorrelation function P{t) fl26p . 

The initial condition (po = also corresponds to a regular trajectory — the periodic 
one-dimensional motion along the x-axis. The autocorrelation function ( |26l) for the cor- 
responding wave packet clearly shows the same periodicity of recurrences, but decreasing 
amplitude of the corresponding peaks reveals much faster spreading than in the former 
case ifo = 7t/2. Clear explanation for it is that any harmonic approximation is no more 
valid for the trajectory, and the Gaussian wave packet is not already a good approxima- 
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tion for the exact time dependent solution. It is worth noting that the wave packet for 
the case (po = undergoes only one-dimensional spreading — along the x-axis, while it 
remains well localized in the y-direction. 

The case with = 7'"/4 sharply differs from the two formers, as it corresponds to a 
chaotic trajectory. Accordingly the autocorrelation function shows almost absent recur- 
rences and even faster spreading of the wave packet. This time the wave packet spreads 
already in two dimensions — over all the chaotic sea. 

Survival probability for the right potential well 



is a natural quantum analogue of the classical quantity N{t)/N{0). The results of nu- 
merical simulations for all the three considered cases are presented on figHH Naturally, 
Pfi = 1 for the trapped state {(fo = 0), as well as for the two others for t < td — 2.5 — 
the classical escape time which is almost the same for both the wave packets. Periodicity 
of the trajectory with ipo = perfectly manifests again in the periodic character of the 
Pnit) for the corresponding wave packet. Graduate decay of the oscillations amplitude is 
due to the spreading. Saturation of the survival probability at value Pr — 0.5 confirms 
that the spreading is one-dimensional in that case: the probability density gets uni- 
formly distributed along the one- dimensional trajectory, which lies between the points 



Xi^2 = ±V^\/l + (recall that E = 1), so exactly one half of l^'p falls into each local 
minimum. 

On the case of the chaotic motion (po = n/A the survival probability is aperiodic and 
quickly saturates at the value pn ~ 0.3. It is in accordance with our arguments of fast 
spreading of the wave packet to uniform two-dimensional distribution of the probability 
density |\E'p over the chaotic sea — indeed, exactly about 30% of the chaotic sea amounts 
for the right minimum, as it can be seen from the Poincare section. If we calculate the 
ratio of the chaotic sea in the right minimum to the whole area of the sea in Poincare 
section, the exact value would be 0.39. 

Our preliminary considerations of three principal types of initial conditions show gen- 
erally good coincidence of the quantum and classical results for the decay of the mixed 
state. Specifically quantum effects in the escape problem, such as the resonance barrier 
penetration and the chaos-assisted dynamical tunneling require further analysis, which 
will be published elsewhere. 

6 Conclusions 

Investigation of the escape from localized areas of configuration space in the existence of 
the mixed state is presented. When the mixed state is present in the system, it is possible 
to "trap" given number of particles in the well. We have considered two possible initial 
distributions. 

For uniform distribution escape law splits into three sections. First section corresponds 
to linear decay, second - to exponential and third forms the plato, which corresponds to 
trapped particles. Number of trapped particles depends only on energy. 
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In the point ensemble case the escape is a two-stage process and number of trapped 
particles depends not only on energy, but on coordinate of injection point too. Only 
trajectories with direction of initial momenta in some cone could escape. Angle of escape 
cone is connected with linear part of regular island in the Poincare section and number 
of trapped trajectories. Moreover, only edge of regular island is necessary to compute the 
number of trapped particles. 
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